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We study the multifractal moments of the current distribution in randomly diluted resistor 
networks near the percolation treshold. When an external current is applied between to terminals 
X and x' of the network, the Ith multifractal moment scales as M^'"'^ {x,x') ~ |a: — x'\^'^'^ , where 1/ 
is the correlation length exponent of the isotropic percolation universality class. By applying our 
concept of master operators [Europhys. Lett. 51, 539 (2000)] we calculate the family of multifractal 
exponents {ipi} for Z > to two-loop order. We find that our result is in good agreement with 
numerical data for three dimensions. 
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I. INTRODUCTION 



Percolation is a leading paradigm for disorder. It provides an intuitively appealing and transparent model of 
the irregular geometry which occurs in disordered systems. Moreover, it is a prototype of a phase transition. Though 
percolation represents the simplest model of a disordered system, it has many applications, e.g., polymerization, 
porous and amorphous materials, thin films, spreading of epidemics etc. 

In particular the transport properties of percolation clusters have gained a vast amount of interest over the last 
decades. Random resistor networks (RRN) are a prominent model for transport on percolation clusters. By means 
of RRN one can study the conductivity of disordered media which might be important for technical applications. 
Nonlinear random resistor networks, for which the voltage drop over an individual resistor depends on some power of 
the current flowing through it, can be exploited to derive various fractal dimensions of percolation clusters. From the 
conceptual point of view, RRN have the advantage that one can formulate a field theoretic Hamiltonian amenable to 
renormalization group analysis. Via RRN one can learn about diffusion on disordered substrates, since the diffusion 
constant D and the conductivity E of the system are related by the Einstein relation 

en 

S = D, (1.1) 

KbT ' ^ ' 

where e and n denote the charge and the density of the mobile particles. The connection of the two problems is 
particularly important, since up to date no direct approach to diffusion on percolation clusters by means of a dynamic 
field theory exists. 

In this paper we study the distribution of currents in RRN. The current distribution has many interesting features, 
one of which is multifractality [^. This means that the distribution is not controlled by one or two relevant length 
scales, but rather by an infinite hierarchy of such length scales. The concept of multifractality was introduced for 
turbulence 1^. It has been applied successfully in diverse areas including diffusion near fractals |4|, electrons in 
disordered media j^] , polymers in disordered media |^ , random ferromagnets |Q , chaotic dissipative systems , and 
heartbeat Q. 

Due to the multifractality infinitely many exponents are needed to characterize the current distribution. Consider 
two connected terminals x and x' of the network. Suppose a current / is inserted at x and withdrawn at x' . The Ith 
moment of the current distribution given (apart from technical details, c.f.. Sec. fl C|) by. 




where the sum runs over all current carrying bonds (the backbone), {■■j^c stands for the average over all diluted 
configurations and ib is an abbreviation for lb/ 1, scales at criticality as 



1 



{x,x') - Ix-x'r"' . (1.3) 

The xi constitute an infinite set of exponents which are not related to each other in a hnear fashion, i.e., the niuhifractal 
moments do not show the usual gap scaling commonly encountered in critical phenomena. 

Each of the m]'-* is associated with a particular subset of backbone bonds having its distinct fractal dimension. 
Let n{i) be the number of bonds carrying current i. Upon applying the saddle point method on finds that the main 
contribution to the ^th moment is given by pO| 

niii)^\x-x'f'^ , (1.4) 

with the multifractal spectrum f{l) and the multifractal exponents xi being related to each other by a Legendre 
transformation. f{l) can be interpreted as the fractal dimension of the subset of bonds dominating Afj'\ 

An elegant approach for studying the multifractal moments is to consider RRN with microscopic noise, i.e., random 
networks in which the conductances of the individual resistors fluctuate about some mean. These noisy RRN were 
originally introduced by Rammal et al. ||ll| to study the effects of Flicker (1//) noise. Flicker noise refers to the low 
frequency spectrum of excess voltage fluctuations measured when a constant current is applied to a resistor. The Zth 
noise cummulant C^^ (cc, x') of the total resistance between the terminals x and x' is proportional to Alf^ {x,x') by 
virtue of Cohn's theorem 112]. 

Historically, the existence of the set of multifractal exponents {xi} was proposed by Rammal et al The authors 
determined several of their exponents for two dimension by numerical simulations. A set of exponents {C2;} equivalent 
to {—xiiy}, where v is the correlation length exponent for percolation, was also proposed by Arcangelis et al. [ p^ . 
Theses authors derived their exponents for several hierarchical structures analytically. The field theoretic description 
of multifractality in RRN was pioneered by Park, Harris and Lubensky (PHL) ||l^. Based on an approach by Stephen 
they formulated a Z?x£'-fold replicated Hamiltonian for noisy RRN. The contributions to the Hamiltonian leading 
to multifractal behavior contain powers of replica space gradients analog to powers of real space gradients, which 



were accounted for as an origin of multifractality by Duplantier and Ludwig |16|. PHL introduced a set of exponents 
{ipn} identical to the set {— x„z^} and calculated it to first order in e = 6 — d, where d denotes the spatial dimension. 
Later on Fourcade and Tremblay [0 gave a reinterpretation of the work by PHL. Batrouni et al. |l8| computed 
several multifractal exponents for d = 3 by numerically solving Kirchhoff's equations. Recently Barthelemy et al. p9| ] 
performed simulations indicating that in the thermodynamic limit the Mj'^ do not exist for / < 0. 

In this article we study the moments of the current distribution by renormalized field theory. We extend our 
real-world interpretation of Feynman diagrams (20|-|22| to RRN with noise. Upon introducing multifractal moments 



for Feynman diagrams we reformulate the field theory of PHL in a way that to our opinion is less complex and more 
intuitive. By carefully analyzing the relevance of the field theoretic operators related to the noise cumulants, we show 
that the multifractality is associated with dangerously irrelevant master operators We calculate the set {ipi} for 
I > to second order in e. Finally, we compare our result to numerical simulations. 

II. THE MODEL 

This section provides background on noisy RRN. It is guided by the work of Stephen jlSj and PHL [ p^ . 

A. Random resistor networks 

Consider a d-dimensional lattice, where bonds between nearest neighboring sites i and j are randomly occupied 
with probability p or empty with probability 1 — p. Each occupied bond has a conductance aij. Unoccupied 
bonds have conductance zero. The bonds obey Ohm's law 

<^^AVJ-Vi)^kj , (2.1) 

where lij is the current flowing through the bond from j to i and Vi is the potential at site i. 

Suppose a current / is injected into a cluster at site x and withdrawn at site x' . The union of all sites belonging to 
all self avoiding paths between x and x' is refered to as the backbone between x and x' . The power dissipated on the 
backbone is by definition 

P = I{V^-V^>) . (2.2) 
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Using Ohm's law, it may be expressed entirely in terms of voltages as 



(2.3) 



Here R{x,x') is the total resistance of the backbone, the sum is taken over all nearest neighbor pairs on the cluster 
and {V} denotes the corresponding set of voltages. As a consequence of the variation principle 



d 



Ip{{v}) 



one obtains Kirchhoff 's law 



(2.4) 



(2.5) 



where h — I {Si,x — Sj^x' ) an d the summations extend over the nearest neighbors of i. 
Alternatively to Eq. ( ^.3| ) the power can by rewritten in terms of the currents as 

P = R{x, x')I^ = E Pbll = P {{h}) , 



(2.6) 



with {7b} denoting the set of currents flowing through the individual bonds, b = {i,j), and pt, = a^^. Obviously 
the cluster may contain closed loops as subnetworks. Suppose there are currents {/('°°p)} circulating independently 
around a complete set of independent closed loops. Then the power is not only a function of I but also of the set of 
loop currents. The potential drop around closed loops is zero. This gives rise to the variation principle 



d 



5/(loop) 



p 



|^|/(loop)| 







(2.7) 



Eq. (2.7) may be used to eliminat e th e loop currents and thus provides us with a method to determine the total 
resistance of the backbone via Eq. (2.6). 

Since the resistance of the backbone depends on the configurations C of the randomly occupied bonds, one introduces 
an average (■ • ■)c over these configurations. It is important to recognize that the resistance between disconnected sites 
is infinite. Therefore one considers only those sites x and x' known to be on the same cluster. Practically this is done 
by introducing the indicator function x{x,x') which, for a given configuration C, is unity if x and x' are connected 
and zero otherwise. Then the Ith moment of the resistance R with respect to the average (• • ■)c subject to x and x' 
being on the same cluster is given by 



{x{x,x')R{x,x')')^ / {X{x,x'))c 



(2.8) 



B. Noise in random resistor networks 



In the following we consider RRN with noise in the sense that the conductances ah of occupied bonds fluctuate about 
some mean. To be specific, the ui, are equally and independently distributed random variables with mean ct and higher 
cumulants A^'-^^. The distribution function / might for example be Gaussian. Nevertheless, our considerations are 
not limited to this particular choice. In order to suppress unphysical negative conductances, the assumption A^'^ <^ ct' 
is made. In general the backbone resistance will depend on the set of conductances of occupied bonds {(Jb\- Its noise 
average will be denoted by 

{R{x, x')}j = / n '^^'-Z ('^f ) (2.9) 

b 

and the corresponding cumulants by 

(2.10) 

A=0 



{R{x,x'yy;^ = ^ln{exp[Ai?(x,x')]} 
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Both kinds of disorder, the random dilution of the lattice and the fluctuation of the bond conductances about their 
mean ct, influence the statistical properties of the backbone resistance. They are reflected by the moments 



M^^\x,x') = {xix.x') {R(x,x'y]^)^/ {x{x,x'))c (2.11) 
and the cumulants 

C« (x, x') (xix, x') {R{x, x')'}^;'')^ I ixix, • (2.12) 



C. Moments of the current distribution 



The noise cumulants characterize the distribution of currents flowing through the network. This section provides 

a relation be twee n the cj^'' and the moments of the current distribution. 

Equation (2^) dcflnes the noise average as an average with respect to the distribution of the bond conductances 
CTfe. Equally well one might express the backbone resistance in terms of the bond resistances and average over the 
distribution of the pi,. Since the ai, are independently and equally distributed, the pb are distributed by the same 
means. Assume that the distribution function of the deviations Spb = pb — p of the resistance of each bond from its 
average p has the form 



{6pb) = -h 
s 



5pb 
s 



and that 



lim Qs (Spb) = S (Spb) 

s— >0 



(2.13) 



(2.14) 



s is a variable with units of resistance which sets the scale of the distribution. With this form of gs, the nth cumulant 
Vn of 6pb tends to zero as s". This follows from the generating function c (As) of the f„: 



exp [c (As)] = {exp {XSpb)}^ ^ I dy h (y) exp (Xsy) = exp | ^ ^ 



(2.15) 



where w„ — c„s" with c„ being constants. In general \^R{x, x'Y^^j^^ depends on the entire set of cumulants {vn}- 
However, in the limit s — > the leading term is proportional to vi as we will see immediately. Consider the generating 



function C (A) of the cumulants {i?(x, x')'}^'^\ 

exp [C (A)] = / n dSpb gs (Spb) exp [XR{x, x')] 

b 

Expansion of the backbone resistance in a power series in the dpb leads to 
exp[C(A)] l[[dyb h{yb)exp 



xRix,x)+xj: j: dpb,---dpb^ 



Vbi--- Vbk 



(2.16) 



(2.17) 



where R^(x^x') is the resistance when 5pb = for every bond b. Equation (2.17) can be rearranged as 



exp [C (A)] — exp 



XR\x,x') + xY, J2 



d^R{x,x') 



k=2 bi,---M 



kl dpb, ■ ■ ■ dpb^ 



Qk 



. dzb^ ■ ■ ■ dzb^ 



J|exp[c(zb) 



i—lb "Pb I 



exp 



XR°ix,x') + J2{XsyciJ2 



1=1 



dR{x,x') 
dpb 



(2.18) 
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where fi are functions of As\ Hence for Z > 2, 



dR{x,x') 



dpb 



l + 0{s) 



(2.19) 



In the hmit s — > the leading term is 



dR{x,x') 



dpb 



b ^ 



21 



(2.20) 



where we have used Cohn's Theorem Eq. (A5). Upon substitution of Eq. ( 2.20| ) into Eq. ( p.l2| ) one finds for the noise 
cumulants 



dj^\x,x')=viMP{x,x') , 



i.e., the noise cumulant C^^^ is proportional to the Ith multifractal moment 



(2.21) 



(2.22) 



of the current distribution. 



D. Generating function 



Our aim is to determine . Hence the task is to solve the set of Kirchhoff's equations ( |2.5| ) and to perform the 
averages over the diluted lattice configurations and the noise. It can be achieved by employing the replica technique 
[p^ . In order to treat the averages (• ■ ■)c and {• • - jy separately, PHL introduced DxE-iold replicated voltages, 



( v^'-'' ■ ■ ■ vt''' 



(2.23) 



Note from the definitions Eq. (2.12) and Eq. (2.1C) that one has to treat the two averages independently in the 

in) Jji] 

calculation of C)j . In contrast, for calculating Af^ it is not necessary to distinguish between the two averages 
because one could also introduce a composite distribution function 



/comp (^) = (l_p)^(^)+p/(^) 

and a single, say D-fold, replication would be sufficient. 

To construct a generating function for the noise cumulants one introduces 

'^A^x) = exp [i A • , 
where A • V x= J2a 'p=i A'-"''^^ vi"''''' and A/O- The corresponding correlation functions 



G 



(x,x'-S^ = (^i^^ix)^_»ix') 



rep 



(2.24) 



(2.25) 



(2.26) 



are defined as 
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1 



g(x,x'-x) = lira ( \ 



3 



X exp 



/3=1 



/3) _ ^(a,/3) 



= 1 (^i) 



// C 



(2.27) 



(2.28) 



with Vx^'' — (Vx^'^\ • ■ • , Vx'^'^'^ ) , and Z is the normalization 



^ ({-^'^ } ' ^) = / n exp ({V^} , {af ^ } , C^) + y E 



(2.29) 



Note that we have introduced an additional power term ^ J2i ■ This is necessary to give the integrals in Eqs. (2.27) 



and (2.29) a well defined meaning. Without this term the integrands depend only on voltage differences and the 
integrals are divergent. Physically the new term corresponds to grounding each lattice site by a capacitor of unit 
capacity. The original situ ation may be restored by taking the limit of vanishing frequency, w — > 0. 

The integrations in Eq. ( ^.27] ) can be carried out by employing the saddle point method. Since the integrations are 
Gaussian the s add le point method is exact in this case. The saddle point equation is identical to the variation prin ciple 
stated in Eq. (2.4). Thus the maximum of the integrand is determined by the solution of Kirchhoff's equations ( |2.5| ) 
and 



g(x,x'S) = (n jexp 



\/3=l 



(2.30) 



The right hand side of Eq. ( ^.30D may be expanded in term of the cumulants defined in Eq. (|2.10[). This gives 



G 



exp 



'^ (-1/2)' ,^ 



(a) {R{x,xi^ 



(2.31) 



where Ki is defined by 



^^(a)-E 

/3=1 



(2.32) 



We learn that the correlation function G can be exploited as a generating function for the noise cumulants via 



{x{x,x'))cC'^' {x,x') 



d 



(-1/2)" 



G 



(x,x';Tj 



(2.33) 



A = 



Note that a4^' = cjj^^ 



E. Field theoretic Hamiltonian 



Since infinite voltage drops between different clusters may occur, it is not guaranteed that Z stays finite, i.e., the 

limit lim£)^o is not well defined. Moreover, A=0 has to b e exc luded properly. Both problems can be handled 
by resorting to a lattice regularization of the integrals in Eqs. (2.27) and ( |2.29 ). One switches to voltage variables 

0= A9 k taking discrete values on a Z?x£'-dimensional torus, i.e., k is chosen to be an _Dxi?-dimensional integer with 
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-M < < M and fc^"-'^) = A;("'^Wod(2M). A9 = 0m /M is the gap between successive voltages and 9m is the 

voltage cutoff. The continuum may be restored by taking 9m — > oo and A9 — > 0. By setting 9m = doM, M = rn? , 
and, respectively, = 9o/m, the two limits can be taken simultaneously via to ^ oo. Since the voltage and current 

variables are conjugated A is affected by the discretization as well: 



A= AA / , AA A9 = n/M , 



(2.34) 



where I is a Dxi?-dimensional integer taking the same values as k- This choice guarantees that the completeness 
and orthogonality relations 



1 



(2M) 



DE 



exp (iX-9)=S 



X,0 mod(2AfAA) 



(2.35a) 



and 



(2M) 



DE 



J2 ( 



i X - 9 



0,0 mod(2A/Ae) 



(2.35b) 



do hold. Equation ( ^.35 ) provides us with a Fourier transform in replica space. In this discrete picture there are 
2M)^^ — 1 independent state variables per lattice site. Upon Fourier transformation one introduces the Potts spins 



$^ (x) = (2M)-^^ A • ^P^ix) =S^» - {2My 

U \ / A tf •: U X 



DE 



(2.36) 



A#0 



subject to the condition i^) = 0- 



Now we revisit Eq. ( 2.27 ). Carrying out the average over the diluted lattice configurations and the noise provides 
us with the weight exp(— iJrcp) of the average (• • •)rop, 



-ffrop = - In ( < cxp 



i 



^ In m <i cxp 

\/3=l 



f C 

X 



(2.37) 



By dropping a constant term Nb ln(l — p), with Nb being the number of bonds in the undiluted lattice, one obtains 



where 



and 



K 



ln<^ 1 



K^) = fE^' 



(2.38) 



(2.39) 



13=1 



1-p 



exp 



ll = l 



(2.40) 
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In the limit of perfect transport, s ^ 0, K (^6^ goes to its local limit K (^6^ — KS^ «, with K being a positive 
constant. The interaction part of the Hamiltonian reduces to 



(2.41) 



This represents nothing more than the (2M)^^ states Potts model which is invariant against all {2M)^^\ permutations 

of the Potts spins $«. 

e 

In the case of imperfect transport this S^^m)'^'^ symmetry is lost. For finite W and A*^") = 0, K (^0^ is an 

exponentially decreasing function in replica space with a decay rate proportional to Then, for large ct, the 

Hamiltonian Hj-^p describes a translationally and rotationally invariant short range interaction of Potts spins in real 

and replica space with an external one site potential h (^6^ . 

Admitting fluctuations of the resistances, A^"^ > 0, results in breaking the rotational 0{DE) replica space symmetry 

of the interaction part of the Hamiltonian. The Fourier transform oi K (^0^, 



K{X] = 



(2M) 



DE 



exp A • 6* Kft: 



(2.42) 



is expediently evaluated by switching back to continuous voltages. 



d exp -i A ■ 6* In < 1 



exp 



E 

.1=1 



(-1/2)' 



where we have dropped a factor (29m 

d 9 exp 



. Taylor expansion of the logarithm yields a series of terms of the form 



-i X - 9 



-aa'f - Y^h [asf Ki (f^ 



1=2 



(2.44) 



where a the bi are constants of order O (s*') . In addition to the expansion of the logarithm we expand in a power 
series in s. 



(2.44) 



d 9 exp 



-i X ■ 9 ~aa 9 



1=2 



(2.45) 



Here the Pi are homogeneous polynomials of order 21 in A which are a sums of terms proportinal to 



i>2 



such that '^j^ili — I. Completing squares in the exponential in Eq. (2.45) gives 

„2 



(2.45) = exp 



exp 



A 
A 



d 9 exp 



1=2 



-aa 



l + Y.{asfpA"9 



1=2 



2a(7 



a^-^^Pi-r = 



(2.46) 



(2.47) 



where we have omitted multiplicative factors decorating the P/. Due to the homogenity of the P;, Eq. ( 2.47 ) can be 
rearranged as 
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(2.47) = exp 



exp 



«2 




A 






[ 






A 




4aCT 


1 




l + 0(s) 



+ (A) + •• 



Pi' X 



(2.48) 



up to multiplicative factors. By keeping only the leading contributions, one finds that K I Xj can be expanded as 



K 



OO 



T + y Wp \ 



-,2p 



(2.49) 



with 



p ^ a and up, A*^'^ /cr^' being expansion coefficients. 
It is known that the terms Wp A are irrelevant in the renormalization group sense for p > 2 (see, e.g. |2C[] ). From 
Sec. [IF can be inferred that the fp, P/ ^A^ are irrelevant as well. However, the terms proportional to Ki ^A^ are 
indispenseble in studying the noise cumulants; they are dangerously irrelevant. Therefore, we restrict the expansion 
ofK(x) to 



if A 



T + wX +^viKi (Tj 

Z=2 



(2.50) 



with w = wi, and vi = ■ Nevertheless the neglected terms will regain some importance later on since they are 
required for the renormalization of the w/. 

The Ki are homogeneous polynomials of order 21. For I > 2 they are possessing a S O (D)^ symmetry. Thus, 
allowing for A'") > results in loosing the rotational O (DE) in favor of the S O {D)^ symmetry. 

It is worth pointing out that vi/w^ ^ A^'^/ct' ^ s', i.e., the condition s ^ translates into vi <^ w' . Consequently 
one has to take the limit vi before the limit w — s- in calculating the exponents associated with the vi. 

We proceed with the usual coarse graining step and replace the Potts spins (x) by order parameter fields ip ^x, 9^ 

which inherit the constraint 12'^ V 9^ =0. We model the corresponding field theoretic Hamiltonian Ti in the 

spirit of Landau as a mesoscopic free energy from local monomials of the order parameter field and its gradients in 
real and replica space. The gradient expansion is justified since the interaction is short ranged in both spaces. Purely 
local terms in replica space have to respect the full S^2M)"'^ Potts symmetry. After these remarks we write down the 
Landau-Ginzburg- Wilson type Hamiltonian 




(2.51) 



where 



D,E 

/v ( A, V« ) = r + A + u; ^ ■ 

a, 13=1 



1=2 0=1 La 



E 



(2.52) 



In Eq. (2.51) we have neglected terms of order Lp"^ or higher which are irrelevant in the ren orma lization group sense. 
T, w and vi are now coarse grained analogues of the original coefficients appearing in Eq. ( 2.50 ). Note again that Ti 
reduces to the usual {2M)^^ states Potts model Hamiltonian by setting vi — Q and w = as one retrieves purely 
geometrical percolation in the limit of vanishing vi and w. 



F. Relevance of the noise terms 



Irrelevant variables that cannot be taken to zero because the quantity one is looking at either vanishes or diverges 
in this limit have been given the name dangerously irrelevant variables by Fisher |25|. Later on this notion was 
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introduced into field theory by Amit and Peliti [g6[ . A characteristic feature of dangerously irrelevant variables is that 
corrections due to them determine the asymptotic behavior of quantities with the above property, so that their effect 
is felt arbitrarily close to the transition ||2^. In this section we show that the vi are dangerously irrelevant. They 
are irrelevant on dimensional grounds, i.e., they are associated with a negative naive dimension. However, we cannot 
simply take the vi to zero by appealing to their irrelevance, because the amplitudes of the noise cumulants vanish in 
this limit. 

In the remainder of this article we focus on vanishing frequency, w = 0. Let P denote the set of parameters 
{r, w, vi}. We introduce a scaling factor b for the voltage variable: 0-^ b . By substitution of (p ^x, 
the Hamiltonian turns into 



n 



V (x, b e) 



\ip' (x, b e) K (a, v») ^' (x, b e) + (x. 



b e 



ip X 



(2.53) 



Renaming the scaled voltage variables 9 ^ b yields 

1 



n 



H>' (x, ^ ) , 



P 



„/N 3 



y: {^re)K (a, w^,) (x, ^ ) + (x, ) 



(2.54) 



Obviously the voltage cutoff is affected by the scaling as well: 0m b0M- However, if the limits are taken in the 
appropriate order, namely D — > and then m — > oo, the dependence of the theory on the cutoff drops out. Thus, we 

can identify and and hence 



H 



where P' — {r, 

Now consider correlation functions 



(x, b ^) 



P 



n 



v{^,e] ,p' 



(2.55) 



Gn [{x,e};T,w,{vi}^ ^ y"p(^(p(xi,0i)---¥.(xAr,0Ar)exp(-H[(^(x,0),p]) , (2.56) 
where Dip indicates an integration over the set of variables |(p ^x, 0^ | for all x and . Equation (2.55) implies 

Gat ({x, 0} ; T, w, - Gn ({x, ; r, b'^w, . (2.57) 



The two-point correlation function G2 is the Fourier transform of ('ijj^{x)'ip "(x) ) . We deduce from Eq. ( 2.31| ) that 

\ A ""A / "hi 

Ki (a) C^;^ {{x,^)-T,w,{vk}) = b-^'Ki (a) Gj^^ ((x,x') ;r,62^, {b^^vu}) . (2.58) 
We are free to choose b^ — w~^. This gives 

Gf ((x,x');r, ^,K}) ((x,x');r,{ ^}) , (2.59) 



where // is a scaling function. We learn from Eq. (2.59) that the coupling constants appear only as Vk/w^. 
Dimensional analysis of the Hamiltonian shows that w \ ^ y? and VkKk ^A^ ^ where /i is an inverse length 

scale, i.e., w A and VkKk i^Xj have a naive dimension 2. Thus Vk/w^ ^ /i^^^*^ and hence the Vk/w^ have a negative 

naive dimension. This leads to the conclusion that the Vk are irrelevant couplings. 

Though irrelevant, one must not set tj; = in calculating the noise exponents. In order to see this we expand the 
scaling function // in Eq. ( ^.59| ), 



G«((x,x');r,.,K}) = -'K^ + G«^ + ...} , 



(2.60) 



with G^'"* being expansion coefficients depending on x, x', and r. It is important to re cogn ize that G^'^j = because 
the corresponding terms are not generated in the perturbation calculation. Equation (2.60) can be rewritten as 
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((x, x') ; r, w, W}) - VI |c« + + • • j , 

(0 



(2.61) 



where the first term on the right hand side gives the leading behavior. Thus C)j vanishes upon setting vi = Q and 
we cannot gain any further information about . In particular we cannot determine the associated noise exponent. 
In other words, the vi are dangerously irrelevant in investigating the critical properties of the Cd"^''. 



III. RENORMALIZATION GROUP ANALYSES 



A. Diagrammatic expansion 



The diagrammatic elements contributing to our rcnormalization group improved perturbation calculation are the 
three point vertex —g and the propagator 



1 -(5» » 

A: 



1 



A,0 



(3.1) 



Note that we have switched to a yp, A j -representation by employing Fourier transformation in real and replica space. 
The notation in Eq. (3.1) is somewhat symbolic. To treat the irrelevant terms viKi properly, we have to expand 



the propagator in a power series in the vi and discard all contributions of higher than linear order in the . In other 
words: the irrelevant terms have to be treated as insertions. ^ 

Eq. ( |3.1[ ) shows that the principal propagator decomposes into a propagator carrying A's (conducting) and one not 

carrying A's (insulating). This allows for a schematic decomposition of principal diagrams into sums of conducting 
diagrams consisting of conducting and insulating propagators. To two-loop order, we obtain the conducting diagrams 
listed in Fig. 1. 



B. Multifractal moments of Feynman diagrams 



From the decomposition in Sec. Ill A a real- world interpretation of the conducting Feynman diagrams emerges 
pO| , pT| . They may be viewed as resistor networks themselves with conducting propagators corresponding to conductors 
and insulating propagators corresponding to open bonds. The parameters s appearing in a Schwinger parametrization 
of the conducting propagators, 



ds exp 



„2 



-S \Tp^ + T + W \ -t- 



1=2 



(3.2) 



correspond to resistances and the replica variables i A to currents. The replica currents are conserved in each vertex 
and we may write for each edge i of a diagram, Ai = Ai ^A, I'*}) j where A is an external current and |k;| denotes a 
complete set of independent loop currents. 

The real- world interpretation suggests an effective way of computing the conducting diagrams. We learn from the 
discussion above, that the irrelevant terms have to be treated by means of insertions 



o^'^ = -\^ij E (^) (p' ^) (-P' - 



A 



(3.3) 



where (p ^p, A^ denotes the Fourier transform of tp ^x, 0^ . The resulting diagrams are of the type displayed on the 
left hand side of Fig. 2. We express the current dependend part of such a diagram in terms of its power P, 



SiVi ^ Ki 



exp 



w 



„2 
J 



-SiVi ^ Ki (^Ai) 



exp 



wP 



(3.4) 
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The summation is carried out by completing the squares in the exponential. The corresponding shift in the loop 
currents is given by the minimum of t he quadratic form P which is determined by a variation principle completely 
analogous to the one stated in Eq. (^). Thus, completing of the squares is equivalent to solving Kirchhoff's equations 
for the diagram. It leads to 



^ind 



s^Vl ^ Ki I Xi +^Cij ({s}) Kj j exp - wR{{s}) A -w'^Bij {{s}) Ki ■ Kj 



(3.5) 



^ind 



Ai — Ci ({s}) A is the current induced by the external current into edge i. Ci {{s}) and Ci_j {{s}) are homogeneous 
functions of the Schwinger parameters of degree zero. Bij {{s}) and the total resistance of the diagram R {{s}) 
are homogeneous functions of the Schwinger parameters of degree one. By a suitable choice of the the matrix 



constituted by the Bij is rendered diagonal, i.e., Bij ^ 
currents and to replace the summation by an integration. 



-'1,3 ■ 



At this stage it is convenient to switch to continuous 



E 

{^} 



(3.6) 



where d k is an abbreviation for Y[a 'p=i dn^"'^') and L stands for the number of independent conducting loops. This 
integration is Gaussian and therefore straightforward. In the limit — > one obtains 



s^viKi ^ + • • • = -s,c, {{s}f viKi (a) + 



(3.7) 



The terms neglected in Eq. (3.7) are not required in calculating the -(/'/■ This issue is discussed in detail in Sec. IIIC. 
Diagrammatically, the calculation scheme can be condensed into Fig. 2. Appendix ^ illustrates the calculations in 
terms of an example. 

So far we have inserted O*-'-* only in one of the conducting propagators. However, each of them has to get an 
insertion. Moreover, the integrations over loop momenta and Schwinger parameters remain to be carried out. All in 
all, each diagram can be written as 



/ (p^ X) - IP {p') - Iw {p') w f -I^^ (p^) viKi (a) + • ■ • 



D{p\{si]) 



(3.8) 



Here D (p^, {si}) stands for the integrand one obtains upon Schwinger parametrization of the corresponding diagram 
in the usual (f)^ theory. C^'-* ({sj}) is defined as 



c« ({..}) = E^^^'(W)'' = E^M A. / A 



.^ind 



21 



(3.9) 



where the sum runs over all conducting propagators of the diagram. Notice the analogy of the C*-'' ({sj}) to the 
generalized multifractal moments we introduce in App. Thus, we refer to the C^'-* ({•Si}) as multifractal moments 
of conducting Feynman diagrams. 



C. Renormalization and scaling 



By employing dimensional regularization and minimal subtraction we proceed with standard techniques of renor- 
malized field theory |^^ . The renormalization of the vi, however, involves some peculiarities that we will discuss in 
this section. 

An operator Oi of a given naive dimension [Oi] inserted one time in a vertex function generates in general new 
primitive divergencies corresponding to all operators of equal or lower naive dimension. Thus, one needs these newly 
generated operators as counterterms in the Hamiltonian. 
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The operators of lower naive dimension can be isolated by additive renormalization, 



Dimensional regularization in conjunction with minimal subtraction leads to Xi^j containing at least a factor r. These 
Xij vanish at the critical point. Hence, the operators of lower naive dimension will not be considered in the following. 
As argued in Sec. Ill B the term proportional to vi in Eq. ( |3.8| ) is generated by inserting the operator O''' . Inserting 
into a diagram with n external legs, see Fig. 3, generates primitive divergencies which must be cancelled by counter 
terms of the structure 



(a) p^-^^ 



P,A , (3.11) 



where 

i 

with iri = r, is a homogeneous polynomial of degree 2r. Note, that the notation we use here and in the following 
is symbolic, since such a counter term has to depend on the entire set of external momenta and currents. ^p, 

is for example an abbreviation for Y\^^i 4> {vi, Xij ■ 

The leading contribution comes from operators having the same naive dimension as O^'-*, i.e., those satisfying 

2 (/ + 2 - 3) = 2 (r + a + 71 - 3) . (3.13) 



Here we expressed the naive dimension with help of Eq. (C8). For n — 2 one is led to ^ > r, i.e., the insertion of 
0(0 generates operators containing homogeneous polynomials in the replica currents of degree equal or lower 21. In 
particular O^'-* generates an operator of type 



viKi {x) 4> {pS)" ■ (3.14) 



The important question n ow is , if the other operators generated by O*^'^ generate operators of this type, too. Consider 
n > 3. With help of Ec^. ( 3.13| ) one obtains ^ — 1 > r > 1, where the second inequality is a consequence of the limit 

Z? — !■ 0. Bearing in mind that max imal homogeneous polynomials of degree Z — 1 in A are generated, we reinsert these 
operators of the type in Eq. ( [3.1l| ) with n > 3 into two-leg diagrams, see Fig. 4. The resulting terms are of the form 



2 



(a)p'">(p,a) , (3.15) 

with the leading contributions satisfying r-\-a + n — i — r' + a' — 1. Thus, r' > r + a — a' + 1, i.e., the homogeneous 
polynomials in A may have a higher degree than 21. However, they are of the type 

Pr' (a) - (a)' n ^' i^T ' (3.16) 

2<i<r 

with iri < r < I — 1 and iri -\- s = r' . These polynomials have a higher symmetry than the original Ki. 

We conclude that O*^'^ generates itself and an entire family of new operators but these in turn do not generate O^'^ 
In principle, the entire family of operators associated with O^'^ has to be taken into account in the renormalization 
proceedure, leading to a renormalization in matrix form 

a^'Ud^'^ =^wa^^ (3.17) 

The vector 

a^'^ = fo('),df (3.18) 
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contains the family associated with 0^''\ For the remaining renormahzations, we employ the same scheme as in [EOl, 



l/2„l/2„e/2 



(3.19a) 
(3.19b) 



In Eq. (3.19) e stands for 6 — d and the factor = (47r)^'^/^r(l + e/2), with T denoting the Gamma function, is 
introduced for convenience. 

According to the arguments given above the rcnormalization matrix 



Z(') = 1 + 0(m) 



(3.20) 



has a particularly simple structure, 



/zw O ••• o\ 
o ••• o 

V o ■■■ o J 



(3.21) 



1 stands for the unit matrix and O symbolizes elements that we do not evaluate. In this paper, we determine Z^^^ to 
the order of two loops. Z, Z^ and Z„ are the usual Potts model Z factors. They have been computed to three-loop 
order by de Alcantara Bonfim et al . Z^^ is known to two- loop order |^ . 

The unrenormalized theory has to be independent of the length scale introduced by renormalization. In 

-(0 

particular, the connected N point correlation functions with an insertion of O must be independent of /x, i.e.. 



(3.22) 



for all N. Eq. ( 3.22 ) translates via the Wilson functions 



/3(«)=A^^ 



K (u) — ^1 



din: 



C(u) = A* 



d In w 



9/1 

7(') (u) 



Jiu)=^^-lnZ 



(3.23a) 
(3.23b) 
(3.23c) 



where the bare quantities are kept fix while taking the derivatives, into the Gell-Mann-Low renormalization group 
equation 



' d d a N ' 

Oj-i ou OT aw 2 



i + 7 



(0 



X G 



N 



X, w A f- ;t, u,/i ) =0 



(3.24) 



The particular form of the Wilson functions can be extracted from the renormalization scheme and the Z factors. At 
the infrared stable fixed point u* , determined by /3 (u*) = 0, the renormalization group equation reduces to 



Ojj, OT ow 2 



1 + 7 



(0» 



X G 



N 



X, w A 



(3.25) 



where 7* = 7 (u*), k* — k (u*), C* — C ("*)i and 7*^')* = 7*^') (u*). 
The matrix 7^'' inherits the simple structure of Z}-^\ 
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( ^{i) o ... o\ 

o ••■ o 



(3.26) 



\ O ■■■ O J 

Owing to this structure, |1) — (1,0, •• ■ ,0)"^ is a right eigenvector of 7*^')* with eigenvalue 7'''*. The remaining right 
eigenvectors with eigenvalues 7^'''*, k > 2, we denote by |fc). The left eigenvectors of 7*^'^* are (1| = (1, O, • • • , O) and 
(fc| = (0, O, • ■ • , O). In terms of the eigenvectors, 7'-'-'* can be spectrally decomposed into 



Now it is important to realize, that 



(0* 
7fe 



(3.27) 



k>2 



k>2 



s(0 



(3.28) 



k>2 



i.e., only A^-^^ contains the operator O^'^ we are interested in. We substitute Eq. (3.27) into the renormalization group 
equation ( |3.2[ ) and act on the entire equation with (1|. The result is 



d ^ d d N 



dw 2 



Gn < X, w a > ;t, u*,^ 



Equation (3.29) is solved by the method of characteristics. The solution reads 



(3.29) 



(3.30) 



AW 



To derive a scaling relation for the correlation functions, a dimensional analysis remains to be performed. It yields 

"'^ ' ^ (3.31) 



Gn [^|x, w a I ;t, u,^ 



■,fl ^T,U,l 



AW 



From Eqs. (3.30) and (3.31) we deduce the scaling behavior 



Gn [ i^^w X } ;t,u,^j. 



AW 



id~2+n)N/2-^j.Q^ gx, g-t'/^w X \ ; g''/-T, u\ti 



(3.32) 



AW 



rj — ^* and 1/ — (2 — k*) ^ are the well known critical exponents for percolation. They are known to third order in e 
P9i: 



1 



206 



'7 = 7= -Trre 



93619 256 



and 



21 9261 



1 5 



2 84 37044 



589 



8168202 7203 



716519 



C(3) 



89 



130691232 7203 



C(3) 



(3.33) 



(3.34) 



Note that C, in Eqs. ( 3.331 ) and (3.34) stands for the Riemann zeta function and should not be confused with the 
Wilson function defined above. 4> = v {2 — C,*) is the resistance exponent known to second order in e [^0|j2^, 



15 



3087 



(3.35) 



The noise exponents V'z are defined by -0; — v (2 



■7 



(0* 



) . The expansion of tpi to second order in e is given below. 
Now we are in the position to derive the scaling behavior of Cj^'^ From Eq. ( 3.32 ) we find upon choosing g = |x— x'|^^ 



and Taylor expanding that the two point correlation function G — G2 scales at criticality as 



G 



(x,x'; a) 



= X — X 



wX |x-x'|^/'^ 



viK, 



(a) |x-x' 



(3.36) 



where we have dropped several arguments for notational simplicity. With Eq. (2.33) the desired scaling behavior of 
cj^^ is now readily obtained as 



C 



(0 



(3.37) 



At this point, we emphasize once more the outstanding role of the O*-''' , which warrants calling them master operators 
p3[ . Each multifractal moment Afj'-* has a master operator as field theoretic counterpart. The master operators are 
highly and dangerously irrelevant in the renormalization group sense. Therefore, each master operator needs in general 
a myriad of other irrelevant operators for renormalization. However, the renormalization of these servant operators 
does not induce their master. It follows, that the servant operators can be neglected in determining the scaling index 
of their master operator, i.e., one is spared the computation and diagonalization of giant renormalization matrices. 

Our e-expansion result for the noise exponents reads 

^ 7 (1 + /) (1 + 21) ^ 12348 (1 + if (1 + 20' 



|313 - 6727 + z|3327 - 40327 - 8;|4 (-389 + 2737) 
I [-2076 + IOO87 + I (-881 + 3367)] I 



672 (1 + if (1 + 2lf *(1 + 2l)\ + (e^ 



(3.38) 



in agreement to first order in e with the one-loop calcu lation by PHL. 7 ~ 0.577215... denotes Euler's constant and ^' 
stands for the Digamma function [Q. Equation (3.38) is valid not only for I > 2 since it can be continued analytically 
down to / = 0. A plot of V'z versus e is gi ven in Fig. 5. We point out that Eq. (3.38) evaluated at Z = 1 is in conformity 
with the result for (j) stated in Eq. ( |3.35 ), i.e., our result for ijji satisfies an important consistency check steaming from 



C 



(1) 



M 



(1) 



Blumenfeld et al. 



proved that V'/ is a convex monotonically decreasing function of I. Note from 



Fig. 5 that our result for captures this feature for reasonable values of e. It reduces to unity in the limit I — *■ 00 



as one expects from the relation of ipoa to the fractal dimension of the red bonds (see App. 
continuation of -0; to ^ = shows that i/'o = vDb up to order O (e"^) as expected (see App. [Gj 



Moreover, analytic 



IV. COMPARISON TO NUMERICAL DATA 



In this section we compare our result for ijji to numerical values. Instead of working with 0; directly we compare 
Tpi/v because data f or ex ponents of this type is available in the literature. For the comparison it is not sufficient to 
simply evaluate Eq. (3.38) at e = 3 or e = 4. The pure e-expansion gives for small spatial dimension poor quantitative 
predictions. However, one can improve the e-expansion by incorporating rigorously known features. We carry out a 
rational approximation which takes into account that iIji/v = 1 in one dimension. Practically this is done by adding 
an appropriate third order term to the e-expansion of liJi/v. We refrain from stating the so obtained formula explictly 
because it is a little lengthy. Instead we plot it for e = 3, i.e. d = 3, in Fig. 6. Our analytic result shows remarkable 
agreement with the the available numerical data for d = 3 js^jl^. For I = our result lays slightly outside the error 
bars of the simulations. However, the deviation of the values is less than 3%. For I = 1,2,3 our result is within the 
error bars of the simulation. There are also numerical values available for d — 2 [^,0. Here, however, the agreement 
is much less pronounced. For ^ = 0, 1, 2, 3, 4 we find a deviation of the order of 30%. It appears that the dependence 
of ipi on dimensionality is too rich in structure to be approximated well at d = 2 by a series of a few terms. 
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V. CONCLUSIONS 



We studied the multifractal moments of the current distribution in RRN by renormalized field theory. Our approach 
thrived on two cornerstones. First, the Feyman digrams for RRN can be interpreted as being resistor networks 
themselves. In this paper we extended our real- world interpretation by introducing multifractal moments for Feynman 
diagrams. The real-world interpretation proves to be a powerful tool which allows a highly efficient calculation of the 
diagrams. 

The second cornerstone was our concept of master operators. Whereas the field theoretic operator associated with 
the resistance exponent (fi is relevant in the renormalization group sense, the operators associated with the ipi>2 are 
dangerously irrelevant master operators. Due to their irrelevance the master operators generate a multitude of other 
irrelevant operators, the servants, which in principle must all be taken into account in the renormalization procedure. 
The servants however, do not influence the scaling index of their master. Without this property, one would have to 
compute and diagonalize entire renormalization matrices for determining the ipi. These renormaliztion matrizes are 
giants for large /. Without the master property it would be practially impossible to compute the "0; for arbitrary I. 

To our knowledge this is the first time that an entire family of multifractal exponents has been calculated to two- 
loop order, at least for percolation. Our result is for dimensions near the upper critical dimension 6 the most accurate 
analytic estimate for the ipi that we know of. It fulfils several consistency checks. Moreover, it agrees remarkably 
well with numerical data for d = 3. As one expects, the agreement suffers by further decreasing the dimension. The 
dependence of the ^pi on dimensionality appears to be too complex for being approximated well at d = 2 by a series 
of a few terms. 

We expect that our concept of master operators can be applied to other systems showing multifractality. It works 
for example to describe the moments of the current distribution in random resistor diode networks. A two- loop 
calculation of the corresponding family of mutifractal exponents will be reported in the near future |3^. Another 
example for the applicability of the concept of master operators is the problem of diffusion near polymers. For this 
problem von Ferber and Holovatch pO| formulated a field theory which comprises dangerously irrelevant operators. 
Due to the symmetry properties of their operators no other irrelevant operators are generated in the perturbation 
calculations of these authors. Thus, the operators studied by von Ferber and Holovatch are particular simple instances 
of master operators. Their scaling index is not influenced by any other operator simply because the number of their 
servants is zero. In this sense these master operators may be called poor. 

It might turn out that the field theoretic operators associated with multifractal quantities are in general master 
operators. In this case the concept of master operators would be a key in understanding the origin of multifractality, 
at least from a field theoretic point of view. 



ACKNOWLEDGMENTS 



We acknowledge support by the Sonderforschungsbereich 237 "Unordnung und grofie Fluktuationen" of the Deutsche 
Forschungsgemeinschaft . 



APPENDIX A: GENERALIZATION OF CORN'S THEOREM AND GENERALIZED MULTIFRACTAL 

MOMENTS 



Consider a generalized power of the form 

b 

where F is some function of bond currents /(,. As argued in Sec. ||, is in general a function of the external current / 
and a complete set of loop currents {^^'^}- We can exploit the variation principle ( pT7| ) to eliminate the loop currents. 
As a result we obtain the as a function of / and {pb} only. The solutions we denote by I™'^. The derivative of the 
so obtained power P {{l™^}) with respect to a bond resistance pb reads 

dpb -^('^ )+4-4-^^'^7(^^- ^^'^ 

The second term on the right hand side vanishes by virtue of the variation principle (I 
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— diU) — -z.^" Old) -0- (A3) 



Renaming 7^ = J™'^ we finally obtain 

dP{{h}) 



F{h) (A4) 



dpb 

as a generalization of Cohn's theorem. For F (Ii,) = one retrieves the original theorem 



dR{x,x') fli 



dpb V I 



2 



(A5) 



Having generalized the power it is natural to generalize the multifractal moments as well. Consider the cumulants 
of the generalized power, 

= ^ln{exp(AF)}^ • (A6) 
In analogy to the resistance cumulants one finds for the leading behavior in the limit s — > 0, 



{Spt} 



/ 



=5:7^(/,)"{5pnf . (A7) 

b 

For networks like those described in Sec. Hone has {Sp'^jf = Vn- However, in a more general situation the individual 
bonds may be composed of a series of elementary resistors. The elementary resistors are assumed to have independently 
and identically distributed resistances with mean p and higher cumulants u/>2- Then 

{SPb}? = ribVn = , (A8) 

^ P 

where Ub denotes the number of elementary resistors constituting bond b and is the average resistance of that bond. 
Upon incorporating a factor 'p~^ into the constants w„ we finally find 

{P''}f =nnY,-p,F{IbT , (A9) 

6 

with X^bPb^ (-^b)" being the nth multifractal moment of F (lb)- 



APPENDIX B: COMPUTATION OF DIAGRAMS I 



In this appendix we illustrate the calculation scheme sketched in Sec. HI B at the instance of diagram A. For the 
sake of simplicity, we focus on its contribution to the renormalization of V2- We neglect all other parts and obtain 



A = 



dsids2-D (p^, {sj) ^exp wP (^,k^ 

K 

^SiV2K2 + S2V2K2 ('^ + ^) } ' 



(Bl) 



where P ^A, = — si k — S2 {^k- + \ j is, according to the real-world interpretation, the power of the diagram and 
D (p^, {si}) stands for 

-,2 



D{p\{s..}) = y / exp 



[si + S2) T - siq- - S2 (q + p)" 



(B2) 
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with being an abbreviation for (27r) ''^^ f d'^q. It is convenient to switch back to continuous currents and replace 
the summation over the loop current by an integration, 

A = — / dsids2D (jpi^ , {si}^ / c? K exp 

Jo J ~oo 

X i S1V2K2 (k) + S2V2K2 + a) > . (B3) 



wP 



The integration over the loop current is simplified by completing the squares in the exponential. One looks for the 
minimum of the quadratic form P I X, k \ . The minimum is determined by a variation principle completely analogous 



to the one stated in Eq. (2.7). Thus completing the squares is equivalent to solving Kirchhoff's equations for the 
diagram. We obtain 



X < S1V2K2 K 



dsids2D (p^, {sj}) exp 

S2 



-R{{s,})w X 



d K, exp 



(Sl ^- S2)w H 



Si + S2 



X + S2V2K2 n 



Sl 



Sl + S2 



(B4) 



where R{{si}) — siS2/(si + S2) is the total resistance of the diagram. Note that S2 X /{si + S2) is, appart from 
a factor i, the replica current induced by the external replica current A into the propagator parametrized by si. 
Sl X /{si + S2) is the replica current induced into the propagator parametrized by S2- In the limit — > we find 



{Airf Jo {si + S2) 



dsids2 



d/2 



exp 



{si + S2)t - R{{s^}) [p^ + w X 



■V2K2(x)+2^1±^^-^X 



(si + S2) 



(.Sl + S2T 



(B5) 



where we have carried out the momentum integration as well. Expanding the exponential and keeping only the terms 
proportional to V2 gives 



A = -9' 



1 



[Anf^ Jo (si + S2) 

Sisi /TA , o^l^i_ ^2 "2 

■ V2^2 A + 2- — — A 



dsids2 



d/2 



exp [- (si + S2) t] 



{si + S2) 

sfsl V2 "2 



(si + S2) 



{S1 + S2T ^ 



— A + w A - 



c2^5 



(Sl + S2) 



V2K2 (a) 



P^ + w A 



(B6) 



The integral over the last term is convergent and therefore neglected. The remaining integrations are rendered 
straightforward by the change of variables si — > tx and S2 — > t (1 — x). Upon expanding the result for small e = 6 — c? 
we obtain 



A = -<, 



e 1 15 \ / 6 w lb w 



p^ + w A 



(B7) 



where we have introduced = (47r)^ '^'^'^ T (1 + e/2) for convenience. We learn, that not only primitive divergencies 
proportional to K2 {^Xj , but also proportional to r A , P^ A and f A j are generated. 



APPENDIX C: SUPERFICIAL DEGREE OF DIVERGENCE OF OPERATOR INSERTIONS 



Consider the insertion of a local operator O by adding a term to the Hamiltonian, 
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H ^ H + J d'^x (x,e^ , 



(CI) 



where O is a local monomial of degree n in the fields ip with A derivatives in real and B derivatives in replica space. 
In a diagram composed of P propagators, V three-leg vertices and the insertion there are 



loops. The topological relation 



L = P-{V + l)-l 



3V + n = 2P + E 



(C2) 



(C3) 



balances the number of legs. Each propagator behaves for large momenta as 1/q^ and hence reduces the superficial 
degree of divergence of the diagram by 2. The insertion increases it by A+B. Thus the superficial degree of divergence 
S [O] of the diagram with insertion is 



5[0]^dL + A + B-2P 



With help of Eqs. (|C2D and (^) one finds 

d\0 



d-6 d-2 2-d 

V + —^n H —E + A + B 



2 2 2 

In contrast, the superficial degree of divergence S of the diagram without insertion is 

, , d-6^^ 
S^d+^V+^E. 

The difference 

S[0]-5= '^-^n + A + B-d 

is identical to the naive dimension [O] of the insertion. For d = 6 it reduces to 

[0] = 2n + A + B-d . 



(C4) 



(C5) 



(C6) 



(C7) 



(C8) 



APPENDIX D: COMPUTATION OF DIAGRAMS II 



Here we give details on the calculation of the conducting Feynman diagrams in listed in Fig. 1. We focus on the 
contributions of the diagrams to the renormalization of the vi, i.e., those terms appearing in Eq. ( ^3.^ ) proportional to 
vi. The other terms appearing in Eq. (3.8) will be omitted throughout the entire appendix for the sake of notational 

simplicity. For details on the calculation of the contributions to the renormalization of w we refer to ||20| . The A- 
independent parts of the conducting diagrams correspond to the usual diagrams found in the literature on the Potts 
model I2II and can be calculated by standard proceedures |28 . 



We start with diagram A. The part of A required in the calculation of ^pi reads for vanishing external momentum 



A = -^viKi ( A / dsids2 / exp [- (si + S2) (t + q^)] <^ si 

^ \ ^ ./n ./r, I \ Si + S2 



where n = 21. Carrying out the momentum integration gives 

^ .00 ^ 



S2 



A = -g^viKi A 



(47r 



.d/2 



dsids2 







(si + S2) 



d/2 



Changing variables, si ^ t (l — x) and S2 — > tx, leads to 



exp [- (si + S2) r] 



S2 



Sl 



Sl + S2 



S1S2 



(si + S2Y' 



(Dl) 



(D2) 



20 



A = -g'^viKi (a) j dx {l~ x) x" ( dt t^-'^''^ exp (-ir) 

^ ^ (47r) ' Jo "'0 



1 



1 



{A^f^ (n+l)(n + 2) 



Expansion for small e — Q — d yields 



r 3- 



G, 



(n+l)(n + 2) e 



(D3) 
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where we have introduced = {Airy'^''^ T (1 + e/2) for convenience. 

The calculation of B is particularly simple. Thus we merely state the result 

B = -g^viKi (a) ^r-/2 . 



(D5) 



Now wc turn to the two-loop diagrams. As an example, we consider the diagram C. As a first step, we determine 
the currents flowing through the conducting propagators. Kirchhoff's law Eq. (2_^) applies to the 4 vertices of the 
diagram. This allows us to eliminate 3 of the 5 unknown currents (one of the vertices is inactive with respect to this 

purpose since the external current A must be conserved). The potential drop around closed loops is zero. Hence we 
can eliminate the two remaining unknown currents and express all currents flowing through conducting propagators 
in terms of the Schwinger parameters and A- The momentum integrations are straightforward. They can be done by 
using the saddle point method which works exactly here since the momentum dependence is purely quadratic. After 
the momentum integration we have 



C = - 



4 / \ 1 
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At this stage, the change of variables si — > ii (1 — x), S2 ^ tix, S3 — » ^2 (1 — y), S4 — > i2y, and S5 — > i3 turns out to 
be useful. It leads to 



C 



—viKi (A T / dtidt2dt3 / dxdy — 

2 V y (47^)" 7o Jo [ti 



exp [-T {h +t2+ h)] 
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The integrations over x and y are straightforward and can be basically looked up in a table 
algebra we obtain 
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After some additional 
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We find it convenient to express the remaining integrals in terms of the parameter integrals given in App. ^. For the 
sake of notational simplicity we introduce the notation 



a—b—c—T 



where /x € {1,3,4}, and 



M^,. (n) = (-1)^+^+'= T&.^M'^ (a, b, c; n) 



da^db^dc'^ 

where /i G {5, 6, 7, 8}. In terms of the parameter integrals we obtain 
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The final result reads 
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4n - 2F2 (n + 3) + 12 2 
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The diagrams D to G can be evaluated in the same fashion. 

As another example we consider diagram H. Determination of the noise cumulants of H leads to 



exp ( -r Y^i^i ■ 
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Here, the change of variables S5 —>■ tix, si — > tiy, S2 ^ ti [1 — x — y), S3 — > ^2, and S4 — > ^3 simplifies the integration. 
We obtain 
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The integrations over x and y are again straightforward. In terms of the parameter integrals we find 



Finally we obtain 
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The diagrams I to L can be treated in a similar manner. 



(D14) 



(D15) 



(D16) 



APPENDIX E: PARAMETER INTEGRALS 

This appendix contains a list of the parameter integrals we use in the calculation of the noise exponents. The results 
stated are obtained by employing the dimensional regularization scheme. The parameter integral given below was 
introduced by Breuer and Jansscn The notation we reserved for a parameter integral we introduced in pc| ] 
but which is not used here. For notational briefness we define 



F,n is related to the Digamma function 'J via 
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In addition to the parameter integrals we use 
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APPENDIX F: CONDUCTING DIAGRAMS IN TERMS OF PARAMETER INTEGRALS 



Here we list our results for the conducting two-loop diagrams in terms of the parameter integrals given in 
For notational simplicity, we show only the parts of the diagrams proportional to viKi 
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APPENDIX G: RELATION TO THE BACKBONE AND THE RED BOND DIMENSION 



From Eq. ( 2.21 ) it is evident that only those bonds with lb — I contribute to C^"*. Consequently, 4'oo is related 
to the fractal dimension d^cd of the singly connected (red) bonds via died = "^oo/v. Conigiio 39 1 proved that 
rfred = 1/j^, which in turn leads to f/^oo = 1. As mentioned above, our result for ipi matches this consistency requirement. 



Another trivial consequence of Eq. (2.21) is, that Cj?"* is proportional to the average number of bonds (the mass) 



of the backbone. Hence ■00 is related to the backbone dimension by 



(Gl) 



This relation can also be obtained on the level of Feynman diagrams. Revisit the definition of the noise cumulants 
for Feynman diagrams Eq. (3.9). In the limit I — > the noise cumulant reduces to the sum of Schwinger parameters 
of conducting propagators, 



(G2) 



Now we take a short detour to our renormalized field theory of diluted networks in which the occupied bond obey a 
generalized Ohm's law V ^ |^,^. In these networks, the nonlinear resistance Rr ix,x') averaged subject to x 
and x' being on the same cluster, 



obeys at criticality 



Mr{x,x') = {x{x,x')Rr{x,x'))(j / {x{x,x'))^ , 
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In analogy to R{{si}) we introduced the notion of the nonlinear resistance i?,- ({si}) of Feynman diagrams. In the 
limit r — » — 1^ we found 

= . (G5) 

i 

Hence we can identify C'-*'-' ({s^}) and ({si}). This leads to the conclusion that ipo = 0-i is related to the 
fractal dimension Db of the backbone via = vDb , and hence we obtain once more Eq. ( pll). Equation (Gl) 
provides another consistency check for our result ( |3.38| ), which is indeed fulfilled. Moreover, Eq. ( |Gi| ) can be used to 
calculate ■00 to third order in e from our three- loop result for Db pl|,p2[, 



We obtain 
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FIG. 1. Conducting diagrams to two-loop order. The bold lines symbolize principle propagators, the light lines stand 
for conducting and the dashed lines for insulating propagators. We point out that the conducting diagrams inherit their 
combinatorial factor from their bold diagram. For example, the diagrams A and B have to be calculated with the same 
combinatorial factor i. 




FIG. 2. Calculation scheme. The hatched blobs symbolize an arbitrary number of closed conducting loops. The solid dots 
indicate insertions. 
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FIG. 3. O''' inserted into a diagram with n external legs. 



Pr(x]p 




-P -A 



FIG. 4. An operator of the type in Eq. (3.11) with n > 3 inserted into a two-leg diagram. 



FIG. 5. Dependence of the noise exponents f/i; on e = 6 — d. 
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FIG. 6. Dependence tpi / v on I in three dimensions. The line shows our rationaly approximated e-expansion. The numerical 
values, symbolized by the dots, are taken from Moukarzel ^] {I = 0) and Batrouni et al. [l — 1, 2, 3). Moukarzel determined 
the backbone dimension Db = tpo/v- Batrouni et al. studied the multifractal moments in a fixed voltage ensemble, i.e., for 
fixed externally applied voltage. The authors state a formula for switching from the multifractal exponents for the fixed voltage 
ensemble to those for the fixed current ensemble. In this formula a minus sign appears to be missing. Correctly, their exponents 
x{n) are related to the iii/v via ipi/v = 2lx{2) — x{2l). 
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